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An enhancement near threshold is observed in the uxj) invariant mass spectrum from the doubly 
OZI suppressed decays of J/i/j — > 7^0, based on a sample of 5.8 x 10 7 J/ip events collected with 
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the BESII detector. A partial wave analysis shows that this enhancement favors J F — + , and 
its mass and width are M = 1812^26 (stat) ± 18 (syst) MeV/c 2 and r = 105 ± 20 (stat) ± 28 
(syst) MeV/c 2 . The product branching fraction is determined to be B(J/ip — » jX) ■ B(X — > u}(f>) = 
(2.61 ± 0.27 (stat) ± 0.65 (syst)) xlO -4 . 

PACS numbers: 12.39.Mk, 13.20.Gd, 13.30.Ce, 14.40.Cs 



QCD predicts a rich spectrum of gg glueballs, qqg hy- 
brids and qqqq four quark states along with the ordinary 
qq mesons in the 1.0 to 2.5 GeV/c 2 mass region. Radia- 
tive J / ip decays provide an excellent laboratory to search 
for these states. Until now, no clear experimental signa- 
tures for glueballs or hybrids have been found. 

Recently, anomalous enhancements near threshold in 
the invariant mass spectra of pp and pK pairs were ob- 
served in J/ip — > 7pp Q] and J/ip — > pKA de- 
cays, respectively, by the BESII experiment. These sur- 
prising experimental observations stimulated many the- 
oretical speculations. Therefore it is of special inter- 
ests to search for possible resonances in other baryon- 
antibaryon, baryon-meson, and meson-meson final states. 

Systems of two vector particles have been intensively 
studied for signatures of gluonic bound states. Pseu- 
doscalar enhancements in pp and lulo final states have 
been seen in radiative J / ip decays jjj 0, HI , and reso- 
nant <p(p structures have also been observed near thresh- 
old in np scattering experiments 0. The radiative J/ip 
decay J/ip — ► jujtp is a doubly OZI suppressed process, 
and its production ratio should be suppressed by at least 
one order of magnitude. Therefore, the measurement of 
this decay and the search for possible resonant states will 
provide useful information on two vector meson systems. 
MARKIII collaboration studied J/ip — > yujcp decays, 
but did not find clear structures in the wcp invariant mass 
spectrum. The final states of ucp were also observed in 
photon-photon collisions by ARGUS kfl llfj experiment 
and the cross sections were measured |l0j . 

In this letter, we report on the measurement of the dou- 
bly OZI suppressed J/ip — > yu<P decay and an enhance- 
ment near threshold in the cu(p invariant mass spectrum, 
using 5.8 x 10 7 J/ip events collected with the upgraded 
Beijing Spectrometer (BESII) at the Beijing Electron- 
Positron Collider (BEPC). BESII is a large solid-angle 
magnetic spectrometer that is described in detail in Ref. 

The J/ip — > joj(p (u> —> tt+tt 7T°, (p —> K+K ) can- 
didate events are required to have four charged tracks, 
each of which is well fitted to a helix that is within the 
polar angle region | cos 6\ < 0.8 in the main drift chamber 
(MDC) and has a transverse momentum larger than 50 
MeV/c. The total charge of the four tracks is required 
to be zero. For each track, the time-of-flight (TOF) and 
specific ionization (dE/dx) measurements in MDC are 
combined to form particle identification Chi-squares for 
the 7r, K, and p hypotheses, and the overall Chi-square 
is determined by adding those of the individual tracks. 



The K+ K~tt + it~ combination is chosen as the combi- 
nation with the smallest combined particle identification 
Chi-square, x 2 ( 7r+7r ~K + K ~), which is required to be 
smaller than x 2 ( 7r+7T - 7r + 7r~) (for the 7r + 7r~7r + 7r~ hy- 
pothesis) and x (K+R-K+K-) (for the K + K~ K + K~ 
hypothesis) to remove the background with 7r + 7r~7r + 7r~ 
and K+K- K+K- final states. 

Candidate photons are required to have an energy de- 
posit in the barrel shower counter (BSC) greater than 40 
MeV, to be isolated from charged tracks by more than 
10°, and to have the difference of angle between the clus- 
ter development direction in the BSC and the photon 
emission direction less than 60°. The number of photons 
is required to be in the range from 3 to 6. 

A five-constraint (5C) energy-momentum conservation 
kinematic fit is made under the J/ip — > ~/K+ K~ir+7{-ir 
hypothesis with the invariant mass of the 77 pair associ- 
ated with the 7T° being constrained to m T o . The combina- 
tion of gammas with the largest probability is chosen as 
the best combination, and events with probability larger 
than 1% are retained. 

To remove backgrounds from J/ip — > K+K ~7r + 7r~7r° 
and J/ip — > K+ K-ir+ir-iT ir , a 5C kinematic fit to 
the J/ip — > K+ K-tt+tt-tt hypothesis and a 6C kine- 
matic fit to J/ip — > K+ K-tt+7T-7t tt (if the number 
of good photons is greater than 4) are performed, and 
the probabilities are required to be less than that from 
the 5C fit to the signal channel. To remove background 
where the 7r° is falsely reconstructed from a high energy 
photon and a second spurious shower, the requirement 
\E 7 i - E 7 2\/\Eji + E l2 \ < 0.90 is applied to the photons 
forming the 7r°. Here, E 7 i and E 7 2 are the energies of 
the two photons. 

Fig. n^a) shows the scatter plot of the m K + K - 
versus the m^,-,! invariant mass after applying the 
above selection criteria. Two clusters are clearly seen, 
which indicate the direct observation of the decays of 
J/ip — y juxp and yepep- The K+K- invariant mass dis- 
tribution is shown in Fig. ^b), where the <p signal can 
be seen clearly. The 7r + 7r~7r° invariant mass distribu- 
tion of candidate events with m^+K- i n th e 4> range 
(\ m K+K- — m <p\ < 15 MeV/c 2 ) is shown as the open his- 
togram in Fig.^c), where to and <j> signals are seen. The 
shaded histogram in Fig. ^c) is the 7r + 7r~7r° invariant 
mass spectrum recoiling against the <p sideband region 
(15 MeV/c 2 < \m K+K - - m | < 30 MeV/c 2 ), where, 
only a very small u> signal is observed which comes from 
backgrounds such as J/ip — > uiK+K~ , uiK*K, etc. Since 
the decays of J/ip —> uxp and ir uj(p are forbidden by C- 
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FIG. 1: (a) The scatter plot of the m K+K - versus the 
7T + 7i- _ 7r° invariant mass, (b) The K + K~ invariant mass 
distribution, (c) The 7r + 7r~7r° invariant mass distribution; 
the open histogram is for candidate events with m K + K - be- 
ing in the </> range, and the shaded histogram is for events 
with m K + K ~ being in the <j> sideband region. (d) The 
K + K~tv + tt~tv° invariant mass distribution for the J/ip — > 
candidate events. The dashed curve indicates the accep- 
tance varying with the u><f> invariant mass, (e) Dalitz plot. 



invariance, the 294 observed LU(f> events present direct ev- 
idence for the radiative J/tp — > jojtfi decay. 

The histogram in Fig.^d) shows the K + K~ir + ir~Tr Q 
invariant mass distribution for events with K + K~ invari- 
ant mass within the nominal <p mass range (\m K + K - — 
to^I < 15 MeV/c 2 ) and the vT + 7r~7r mass within the uj 
mass range (1777^+^-^0 — m u \ < 30 MeV/c 2 ), and a struc- 
ture peaked near oj<p threshold is observed. There is no 
evidence of an n c signal in the uj(p invariant mass spec- 
trum. The dashed curve in the figure indicates how the 
acceptance varies with invariant mass. The acceptance 
decreases as the invariant mass of uxp becomes smaller 
due to the decay of the kaon. The peak is also evident as 
a diagonal band along the upper right-hand edge of the 
Dalitz plot in Fig. ^e) . There is also a horizontal band 
near rrii K+K - = 2 (GeV/c 2 ) 2 in the Dalitz plot, which 
mainly comes from background due to J /ip — > ujK*K. 

To ensure that the structure at the oj(p mass thresh- 
old is not due to background, we have studied potential 



background sources using both data and Monte Carlo 
(MC) data. Non-w and non-0 background are studied 
using u> and <f> sideband events. Fig. [2Ja) shows the 
K + K~ir + TT~ir° invariant mass of events within the uj 



sideband (50 MeV/c 2 < |m„ 



< 80 MeV/c 2 , 



\ m K+K- - m <p\ < 15 MeV/c 2 ), Fig. |2fb) shows the 
corresponding spectrum of events within the <f> side- 
band (Im^-^o - m u \ < 30 MeV/c 2 , 15 MeV/c 2 < 
\™>K+K- - m 4>\ < 30 MeV/c 2 ), and Fig. Efc) shows 
the events in the corner region, which is defined as (50 



MeV/c 2 < K+^o 



< 80 MeV/c 2 , 15 MeV/c 2 < 



\m K + K - — m^l < 30 MeV/c 2 ). The background is es- 
timated by summing up the normalized backgrounds in 
Fig. Eta) and Fig. E{b) and subtracting that in Fig. 
Etc), and it is shown as the shaded histogram in Fig. 
Eld). No evidence of an enhancement near uxp threshold 
is observed from the non-w and non-c/3 background events. 
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FIG. 2: The K + K~ ty + tv~ ty° invariant mass distribution for 
(a) the events in the uj sideband; (b) the events in the <f> 
sideband; (c) the events in the corner region; (d) for events in 
the uj(f) range, as described in the text. The shaded histogram 
in (d) represents the background distribution obtained from 
the sideband evaluation. 

Exclusive MC samples of J/ip decays which have sim- 
ilar final states are generated to check whether a peak 
near u>4> mass threshold can be produced. The main back- 
grounds come from J/tp — > ojK*K,K* — > Ktt°. About 
45+/-17 J/tp — > ujK*K,K* — > Ktt° events remain in 
the u)(j) invariant mass. However, they peak at the high 
mass region and do not produce a peak near the thresh- 
old. We also checked possible backgrounds with a 60 
million Monte-Carlo J/tp — > anything sample, generated 
by the LUND-Charm model None of the MC chan- 
nels produces a peak near threshold in the uj(f> invariant 
mass spectrum. In addition, the data taken at the e + e~ 
center of mass energy of 3.07 GeV, with a luminosity 
of 2272.8 ± 36.4 nb^ 1 are used to check the continuum 
contribution. No events are survived. As a check, the 
measurement of the branching fraction of J/ip —> "f<p4> — > 
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jtt + tt~tt°K + K~ is performed, and the result is consis- 
tent with that from J/tp — > jcpcp — > r )K + K~ K + K~ , but 
with a larger error. 

A partial wave analysis is used to study the spin-parity 
of the enhancement, denoted as X. The amplitudes are 
constructed with the covariant helicity coupling ampli- 
tude method [lflj . and the maximum likelihood method 
is utilized. The decay process is described with sequen- 
tial 2-body or 3-body decays: J/tp — > jX, X — > uxp, 
u> — > 7r + 7T _ 7r°, and — > The resonance X is pa- 

rameterized by a Breit-Wigner with constant width, and 
the background is approximated by non-interfering phase 
space. The uj decay amplitude is not considered in the 
fit. The details of the likelihood function construction 
can be seen in Ref. 0] . 

When J/tp — > jX, X — > uxp is fitted with both the 
oj(j) and jX systems being 5-wave, which corresponds 
to a X = ++ scalar state considering the C-parity 
of the uxp system, the fit gives the best log likelihood 
value of -79.55. The fit gives 95 ± 10 events with mass 
M = 1812126 MeV/c 2 , width T = 105 ± 20 MeV/c 2 , 
and a statistical significance larger than 10 a. Fig. CJa) 
shows the comparison of the uxp invariant mass distri- 
butions between data and MC projection with the fit- 
ted parameters. The comparisons of the angular distri- 
butions between data and MC projection for the events 
with the invariant mass less than 2.0 GeV/c 2 are shown 
in Fig. Eft) - (f). 

If the decay of the 7X system in J/tp — > jX, X — > uxp 
is treated as X>-wave or a combination of both T> and 
<S-wave, the mass and width of the X, as well as the 
log likelihood value do not change much. However, if 
the decay of the X to uxp is fitted with a "D-w&ve, the 
log likelihood value gets worse by about 40. This means 
that the orbital angular momentum of the X — > uxp de- 
cay can be well separated between 5-wave and 2?-wave, 
while it is difficult to determine in the 7A system. A 
fit with P-wave decays in both the uxp system and 7 A 
systems, corresponding to X being a h pseudoscalar 
state, makes the log likelihood value worse by 58. Theo- 
retically, the h hypothesis can be well separated from 
the ++ hypothesis by the distributions of the polar an- 
gle of the K + (9k+ ) in the <p> rest system, the polar angle 
of the normal to the u> decay plan {6 J) in the u> rest sys- 
tem, and x, the angle between the azimuthal angles of 
the normal to the u> decay plane and the momentum of 
a K from <j> decay in the X rest system. We also tried 
to fit the resonance X with 2 ++ and 2 h spin-parity hy- 
potheses with all possible combinations of orbital angular 
momenta in the uxp and 7A systems. The log likelihood 
values of the best fits are -Log(£) = -74.80 and -Log(£) 
= -63.81 for 2 ++ and 2 h assignments, respectively. Al- 
though the fit is not too much worse for the 2 ++ case, 
there is only one free parameter for ++ , while there are 
four free parameters for 2 ++ . Also for 2 ++ , D-wave is 
required in X — » uxp; if only 5-wave is used in the fit, 
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FIG. 3: Comparison between data and MC projections using 
the fitted parameters for the uxp invariant mass distribution 
and the angular distributions for the events with uxp invariant 
mass less than 2.0 GeV/c 2 . Points with error bars are data, 
the solid histogram is the MC projection, and the dashed line 
is the background contribution, (a) The uxj> invariant mass 
distribution; (b) The polar angle of radiative photon (0 7 ); (c) 
The polar angle of the <p in the uxp rest system (#0). (d) The 
polar angle of kaon in the <p rest system (0k) (e) The polar 
angle of the normal to the u) decay plane in the ui system. 
(f) The x distribution - the angle between azimuthal angles 
of the normal to the w decay plane and the momentum of a 
kaon from <p decay in the X rest system. 



-Log(jC) = -67, which is much worse than the final fit. 
Therefore we conclude that the J PC of the enhancement 
X favors 0++. 

Using the selection efficiency of 1.44%, determined 
from Monte-Carlo simulation, we obtain the product of 
the branching fractions as: 

B(J/ip -> 7A) • B(X -> lo<P) = (2.61 ± 0.27) x 10~ 4 . 

Since phase space J/ip — > "/uxp decays exist, fitting 
with an interfering phase space (0 + ) is also performed, 
and the differences between fitting with non-interfering 
phase space for the mass, width, and branching ratio are 
0.5%, 21.9% and 13.8%, respectively. The differences will 
be included as systematic errors. 

The systematic uncertainties on the mass and width 
come from the uncertainties in the background, the mass 
calibration, and the interference with phase space, as 
well as possible biases due to the fitting procedure. The 
latter are estimated from differences between the input 
and output masses and widths from MC samples, which 
are generated as J/tp —* jQ ++ , ++ — > uxp with S- wave 
in both the uxp and jX systems. The uncertainties in 
the background include the uncertainty in the amount of 
background as well as the treatment of the background 
in the fitting. We also tried to subtract the background 
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determined from the sidebands in the fit instead of us- 
ing the non-interfering or interfering phase space back- 
ground, and the differences are taken as systematic er- 
rors. The total systematic errors on the mass and width 
are determined to be 18 MeV/c 2 and 28 MeV/c 2 , re- 
spectively. The systematic errors in the branching frac- 
tion measurement mainly come from the efficiency dif- 
ferences between the Monte-Carlo simulation and data, 
which include the systematic uncertainties of the track- 
ing efficiency, the photon detection efficiency, the particle 
identification efficiency, the kinematic fit, and the to and 
4> decay branching fractions, the amount of background, 
MC statistics, the fitting procedures, different treatment 
of background, and the total number of J/ip events. The 
total relative systematic error on the product branching 
fraction is 25%. 

In summary, the doubly OZI suppressed decay of 
J/ip — > ^u)(f),u) — > 7r + 7T~7r°,0 — > K + K~ is studied. 
An enhancement near uj(j) threshold is observed with a 
statistical significance of more than IOct. From a par- 
tial wave analysis with covariant helicity coupling am- 
plitudes, the spin-parity of the X — ++ with an S- 
wave uj(j) system is favored. The mass and width of 
the enhancement are determined to be M = 1812^26 
(stat) ± 18 (syst) MeV/c 2 and r = 105 ± 20 (stat) ± 
28 (syst) MeV/c 2 , and the product branching fraction is 
B(J/ip -> jX) ■ B(X -> u<t>) = (2.61 ± 0.27 (stat) ± 0.65 
(syst)) x 10~ 4 . The mass and width of this state are not 
compatible with any known scalars listed in the Particle 
Data Gro up (PPG) Il5l . It could be an unconventional 
state [lfiL fl7l TisL ll3.T2(T|. However, more statistics and 
further studies are needed to clarify this. 
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